ENHANCED LIFESPAN OF SMOOTH SOLUTIONS OF A 
BURGERS-HILBERT EQUATION* 

JOHN K. HUNTER t AND MIHAELA IFRIM * 



Abstract. Wc consider an initial value problem for a quadratically nonlinear inviscid Burgers- 
Hilbert equation that models the motion of vorticity discontinuities. We use a normal form trans- 
formation, which is implemented by means of a near-identity coordinate change of the independent 
spatial variable, to prove the existence of small, smooth solutions over cubically nonlinear time-scales. 
For vorticity discontinuities, this result means that there is a cubically nonlinear time-scale before 
the onset of filamentation. 
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1. Introduction. We consider the following initial value problem for an inviscid 
Burgers-Hilbert equation for u(t, x; e): 

u t + euu x =H[u], 
u(0, x; e) = uq{x). 



In (1.1 1, H is the spatial Hilbcrt transform, e is a small parameter, and Uq is given 
smooth initial data. This Burgers-Hilbert equation is a model equation for nonlinear 
waves with constant frequency [T] , and it provides an effective equation for the motion 
of a vorticity discontinuity in a two-dimensional flow of an inviscid, incompressible 



fluid [Il|6j. Moreover, as shown in JT], even though (1.1) is quadratically nonlinear 
it provides a formal asymptotic approximation for the small-amplitude motion of a 
planar vorticity discontinuity located at y = eu(t, x; e) over cubically nonlinear time- 
scales. 

We assume for simplicity that x G R, in which case the Hilbert transform is given 

by 



U[u](t,x;e) = p.v.i J 



u(t,y;e) 
dy. 



We will show that smooth solutions of (1.1) exist for times of the order e as e — > 



Explicitly, if iP(R) denotes the standard Sobolev space of functions with s weak 
L 2 -derivatives, we prove the following result: 

Theorem 1.1. Suppose that uq G i? 2 (R). There are constants k > and eo > 0, 
depending only on \\uq\\h2 : such that for every e with |e| < £o there exists a solution 

ueC (r; H 2 (R)) n C 1 (7 e ; H 1 (R)) 



of (1.1) defined on the time-interval I e = [— fc/e 2 , /c/e 2 ] . 

The cubically nonlinear 0(e~ 2 ) lifespan of smooth solutions for the Burgers- 
Hilbert equation is longer than the quadratically nonlinear 0(e _1 ) lifespan for the 
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inviscid Burgers equation Ut + tuu x = 0. The explanation of this enhanced lifespan 



is that the quadratically nonlinear term of the order e in (1.1 1 is nonresonant for 
the linearized equation. To see this, note that the solution of the linearized equation 
Ut = H[u] is given by u = e* H uo, or 

u(t, x) = uq(x) cost + ho(x) sint, ho = H[uo], 

as may be verified by use of the identity H 2 = — I. This solution oscillates with 
frequency one between the initial data and its Hilbert transform, and the effect of 
the nonlinear forcing term euu x on the linearized equation averages to zero because 
it contains no Fourier component in time whose frequency is equal to one. 

Alternatively, one can view the averaging of the nonlinearity as a consequence of 
the fact that the nonlinear steepening of the profile in one phase of the oscillation 
is canceled by its expansion in the other phase. This phenomenon is illustrated by 
numerical results from Q] , which are reproduced in Figure |1.1| The transition from 
an 0(e _1 ) lifespan for large e to an 0(e" 2 ) lifespan for small e is remarkably rapid: 
once a singularity fails to form over the first oscillation in time, a smooth solution 
typically persists over many oscillations. 




2 It £ 



Fig. 1.1. Logarithm of the singularity formation time T s for the Burgers- Hilbert equation 
versus the logarithm of 2ire for fixed initial data uq. Numerical solutions are shown by diamonds. 
The steeper line is a formal asymptotic prediction from [T^ for £ C 1, which gives T s = 2.37 e~ 2 . 
The shallower line is the singularity formation time for the inviscid Burgers equation, which gives 
T s = t' 1 . (See US for further details.) 

In the context of the motion of a vorticity discontinuity, the formation of a sin- 
gularity in a solution of ( |1.1| corresponds to the filamentation of the discontinuity 
[3]. The result proved here corresponds to an enhanced lifespan before nonlinear 
'breaking' of the discontinuity leads to the formation of a filament. 

There are three main difficulties in the proof of Theorem [hi] The first is that the 



presence of a quadratically nonlinear term in ( 1.1 1 means that straightforward energy 
estimates prove the existence of smooth solutions only on time-scales of the order 
e _1 . Following the idea introduced by Shatah [7j in the context of PDEs, and used 



BURGERS-HILBERT EQUATION 



3 



subsequently by other authors, we remove the quadratically nonlinear term of the 
order e by a normal form or near-identity transformation, replacing it by a cubically 
nonlinear term of the order e 2 . The second difficulty is that a standard normal form 
transformation of the dependent variable, of the type used by Shatah, leads to a loss 
of spatial derivatives because we are using a lower-order linear term H[w] to eliminate 
a higher-order nonlinear term euu x . The third difficulty is that (1.1 ) is nondispersive 



and solutions of the linearized equation oscillate but do not decay in time. Thus, we 
cannot use any dispersive decay in time to control the loss of spatial derivatives. 

The key idea in this paper that avoids these difficulties is to make a transformation 
of the independent variable, rather than the dependent variable. We write 



h(t, x; e) = H[u](£, x; e) 



(1.2) 



and define 



g{t,£;e) =h(t,x;e), x = £ - eg(t, £; e). 



(1.3) 



Then, as we will show, the transformed function g(t, £; e) satisfies an integro-differential 
equation of the form 



gt{t,£;e) = p.v.- / — 



-e J (£ - Qgfat; e)(j> I — ~ ;el d£ 



(1.4) 



where </>(c; e) is a smooth function, given in Lemma 2.1 The term of the order e has 



been removed from (1.4), and the equation has good energy estimates that imply the 



enhanced lifespan of smooth solutions. 



The interpretation of the transformation ( 1.3 1 is not entirely clear. On taking the 



Hilbcrt transform of ( 1.1 1 we get h t = —u + 0(e), so that 



x t 



-m 



-eh t + 0(e 2 



0(e 2 



Thus the transformation £ i— > x in ( 1.3 ) agrees up to the order e with a transformation 



from characteristic to spatial coordinates for ( 1.1 1. The coordinate £, however, differs 
from x even when t — 0, and the use of characteristic coordinates does not appear 
to simplify the analysis. As a partial motivation for (1.3 1, we show in Section [3] that 
it agrees to leading order in e with a normal form transformation of the dependent 
variable that is given in pQ. We were not able, however, to use the latter normal 



form transformation to prove Theorem 1.1 because of the loss of derivatives in the 
higher-order terms. 

We consider (1.1 1 on the real line for simplicity. Equation ( |1.1[ ) is nondispersive 
and solutions of the linearized equation oscillate in time. Thus, our proof does not 
depend on any dispersive decay of the solutions in time, and a similar result would 
apply to spatially periodic solutions. Theorem |1.1| is also presumably true in H s 
for any s > 3/2; we consider s — 2 to avoid complications associated with the use 
of fractional derivatives. A proof of singularity formation for (1.1 1 under certain 
conditions on uq and e is given in [2]. 
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2. Proof of the Theorem. In this sectio n, w e prove Theorem [TTJ It follows 
from standard energy arguments (e.g. [5]) that (1.1) has a unique local 7f 2 -solution 
in a time-interval J £ depending on the if 2 -norm of the initial data and e. Moreover, 
for any s > 2, the solution remains in H s if the initial data is in H 8 and depends 
continuously on the initial data in C(J £ ; H s ). Thus, in order to prove Theorem it 
is sufficient to prove an a priori if 2 -bound for smooth solutions u € C°° (I e ; H°°(M.)) 
where H°°(M.) — D'^ 1 H S (M.). To derive this bound, we first transform the equation 
to remove the order e term and then carry out ff 2 -estimates on the transformed 
equation. The required computations, such as integrations by parts, are justified for 
these smooth solutions. 

2.1. Near-identity transformation. Let h denote the Hilbert transform of u, 



as in (1.2). Taking the Hilbert transform of (1.1), using the identity 



H [u 2 - h 2 } = 2hu 
and the fact that u = — H[/i], we find that h satisfies the equation 
h t + e{H [hh x ] - hH [h x ] - H [h] h x } = H[h}. 



(2.1) 



We will make the change of variables 41.31) in (2.1), so first we discuss (1.3). 



The map £ i— > x is smoothly invertible if \eg^ < 1, which holds by Sobolev em- 
bedding if ||e<?||#2 is sufficiently small. Specifically, we have the Gagliardo-Nirenberg- 
Moser inequality 



llflelk- <N}\g\\% 4 \\g^\ 



1 3/4 
L 2 ' 



where we can take, for example, 



N = 



We assume throughout this section that 



^NH^ 4 ||e^|| 3 L / 4 <-, 



which ensures that 



(2.2) 



(2.3) 



(2.4) 



By the chain rule, 



h,_ r = 



Thus, if (2.4) holds, then 



Hence, since H is an isometry on H s , 



I"IIl 2 = \\9\\li 



5/2 



9tt 



/ h 2 dx = [ g 2 (1 - e 5e ) <%= [ g 2 d^ f h 2 xx dx = f 
JR JR Jr Jr Jr 



(2.5) 



9k 

(1 - eg 6 f 



lf«IU 2 < \\u xx \\l 2 < 2 5/2 ||3alU 2 



(2.6) 
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and ff 2 -estimates for g imply ff 2 -estimates for u. 

Conversely, one can use the contraction mapping theorem on Co (R) , the space of 
continuous functions that decay to zero at infinity, to show that if ho G Cq (R) and 

lleMU- < 1 (2-7) 

then there exists a function go (•; e) € Co(IR) such that 

ho (£ - eg (£;e)) =g (£;e). 

The function go is smooth if ho is smooth, and ||e<?o£lU°° < 1/2 if Heft-ozlU 00 < 1/3. 
Thus, we can obtain initial data for g from the initial data for h. 
From (1.3), we have 



fJi 



l - eg? 

where we use the notation 



H[h] = p.v.- 

7T 



1 



t-Z-e(g-g) 



9 = 9{t,£;e), g = g(t,£;e). 



Using these expressions, together with ( |2.5| , in (2.1) and simplifying the result, we 
find that <?(£,£; e) satisfies the following nonlinear integro-differential equation: 



1 



gt = p-v.- 



g + e(g - 2g)g g - - e 2 (g - g)g^g^ 
S-£-e(g-g) 



d£. 



Subtracting off the leading order term in e from the integrand, we may write this 
equation as 



g t = H[ff] 

where 



7T 



.9 



g-g 



9t 



de (2. 



H[ 5 ](t,£;e)=p. 



1 



gM^de 



denotes the Hilbert transform of g with respect to £ and 

x = eg(t,£,;e), x = £ - eg(t,£;e). 



The integral of the order e in ( 2.8 1 is not a principal value integral since the integrand 
is a smooth function of (£,£)■ Finally, we observe that this equation can be put in 
the form (1.4). 

Lemma 2.1. An equivalent form of equation \27oy is given by 



1 



g t = H [g] - -e*d t jj^ - { |— |; e J d£, 



(2.9) 



whe 



(c;e) = -^{log(l-ec)+ec}. 



(2.10) 



G 
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Proof. First, we check that (2.9) is well-defined. Abusing notation slightly, we 
write 



9-9 



From (2.10), 



^c(c; e) 



1-ec 



(2.11) 



(2.12) 



so \4>{c\ e)| < c when |ec| < 1/2, which is implied by (2.4). In that case 



(£-£)s 5 ~0(c;e) df 



< 



(.9 - fif) % c 



We use |# — #| < 2||<7||£=° in the right hand side of this inequality and apply the 
Cauchy-Schwartz inequality to get 



<2\\g\\ L ^\\ L ,\\c\ 



L- 



(2.13) 



where 



~ N 2 
9-9 

i-l 



1/2 



denotes the L 2 -norm of c with respect to £, which is a function of £. Temporarily 
suppressing the (t; e)-variables and denoting the derivative of g with respect to £ by 
ff'CO = 9i(0i we have from the Taylor integral formula that 

c= y </(f + r(£-£)) dr, 
and the Cauchy-Schwartz inequality implies that 



irds 



\4h (0= / c 2 d| 

| jf j" ({ + r({ - 0) <>' (f + s (f - 0) 

i ' i ' (/, 9 " + - { ') *) ^ (/, ^ ( £ " + * - ^ *) ** 



< 



< 



< 4 



Jo 



drd, 1 



5 ,2 (0de 



Thus, 



1/2 

supl / c 2 d£) <2|| ff5 || i2 



(2.14) 
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Using this estimate in (2.13), we get 



sup 



<%ll 



Thus, the ^-integral in (2.9) converges when <? e iJ 1 (E) and is, in fact, a uniformly 
bounded function of £. 

To verify that (2.9) agrees with (2.8), we take the ^-derivative under the integral 



in (2.9), use (2.12) which implies that 

<t>c (c; e) 



9~9 



and integrate by parts in the result. This gives 



9t = H [g] + -e 2 

7T 



9-9 



x — x 



gc^ - (£ - Og|c c d|. 



Using the equations 



9i 



9e 



(2.15) 



(2.16) 



in (2.15) and comparing the result with (2.8) proves the lemma. □ 



2.2. Energy Estimates. Multiplying (2.9) by g, integrating the result with 
respect to £, and integrating by parts with respect to £, we find that the right-hand 
side vanishes by skew-symmetry in (£,£) so that 



dt 



\9\\& 



0. 



(2.17) 



The conservation of ||<?||z,2 is consistent with the conservation of from ( |1.1[ ). 

Hence, from ( |2~6| ), we have \\g\\ L * = \\9o\\l* = \\uq\\l^- 

Differentiating (2.9) twice with respect to £, multiplying the result by g^, inte- 
grating with respect to £, and integrating by parts with respect to £, we get 



d 
dt 



[ gl e d£ = e 2 I 

Jr 



where 



I = 



(£-0^(qe) 



d£d£. 



(2.18) 



(2.19) 



The following lemma estimates / in terms of the i? 2 -norm of g. 



Lemma 2.2. Suppose that I is given by (2.19) where <fi is defined in (2.10), and 

(2.20) 



c is defined in (2.11). There exists a numerical constant A > such that 

\I\<A\\ge\\ L A9x\\h 



whenever g € H°°(R) satisfies 

Proof. We first convert the ^-derivative in the expression (2.19) for / to a £- 
derivative. Let 



&(c;e)=cf>(c; e), 
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where a prime on $ and related functions denotes a derivative with respect to c. It 
follows from (2.161 that 

(£ - £ ) = (£ - [ c - - £> c i] <^ c ; £ ) 

= (^-^( C ;e)-(^-0 2 $ f -(c; e ). 



We use this equation in (2.191 and integrate by parts with respect to £ in the term 
involving $. Since g is independent of £, this gives 

1= I 9m9l k-e>(c;e)l d£df 



(2.21) 



where 



*(c;e) =c^(c;e)-2$(c;e). 



Expanding the derivatives with respect to £ in (2.21 1, using (2.11 1 to express 
in terms of g^, and integrating by parts with respect to £ in the result to remove the 
third-order derivative of g, we find that / can be expressed as 



where 



h = 
h = 
h-- 

The functions ^>" , ty'" are given explicitly by 



I = --h + 3I 2 -I 3 
vj/" ( c; e )c £ gl^dl 

(e-e>'"(c;e) c ^ K d^e 



(2.22) 



*"(c;e) 



(l-ec)* 



*"'(c;e) 



1 + ec 
(1 -ec) 3 ' 



In particular, if |ec| < 1/2, which is the case if g satisfies (2.3), then 
|*"(c;e)|<4|c|, |*"'(c; e)\ < 12. 



(2.23) 



We will estimate the terms in (2.22) separately. 

Estimating I\: Using (2.23) in (2.22), we get that 

|/ii <4 |cc c 4| de« 

|cc«l ^1 



< 4 sup 

< 4 sup 



c 2 ^ 



1/2 



(2.24) 



1/2" 
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By a similar argument to the proof of (2.14), using Taylor's theorem with integral 
: Cauchy-Schwartz 

9 - 9 - (£ - £)9t 



remainder and the Cauchy-Schwartz inequality, we have from (2.11) and (2.16) that 

2 



10 Jo 

f 1 /•! 



(1 - r)(l - s)g' (£ + r(t - !)) g' (i + s(t[ - £)) dfdrcfe 



< 



(l-r)(l-s) 



< 



o JO 
16 „ 

< yllsdl 



</ 2 (l+r(e-0) < 
' " 1 (l-r)(l-.s) 



' 2 (i+^-o) « 



1/2 



drds 



2 

L 2 ' 



Thus, 



SUP ( / c| G?£ 



1/2 



< 3 ll.^ee II i- 2 - 



Using (2.14) and (2.25) in (2.24), we get that 



\h\ < Ai\\9t\\ L A9ii\\l^ 



where A\ — 32/3 is a numerical constant. 



Estimating I 2 : Using (2.23) and (2.25) in (2.22), we get that 



\h\ <4 / \cclg 66 \ d£dt 



< 4y ( sup r 

64 



< 



\\9xWh I (bupm) | fl «ide. 



Suppressing the (t; e)-variables, we observe from (2.11) that 

9~9 



sup \ c\ = sup 



1 



— SU P 



i(z) dz 



where 



= sup 



le-ei 



Iff'WI^ 



(2.25) 



(2.26) 
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is the maximal function of g' = g^ , denned using intervals whose left or right endpoint 

Using this inequality and the Cauchy-Schwartz inequality in (2.26), we find that 

64„ 



/a < 



Y\m\\L 2 \\9ii\\L^ 



The maximal operator is bounded on L 2 , so there exists a numerical constant M such 
that 



hlh* <M\\gt\\&. 
For example, from [3], we can take 

M = 1 + V2. 

It follows that 

\h\ < M9t\\L49tt\\h 
where A 2 = MM/9. 



Estimating I 3 : Using (2.161 in (2.22), we we can rewrite I 3 as 
Splitting this integral into two terms, we get I 3 = I 3 — I 3 where 



l' 3 = &"(c;e)cc 2 g ee dtd(, I" = I *"'(c; e)cf«^ 



Using (2.23), we have 

\I' 3 \ < 12/ \cc 2 m \dtdi \m<12 [ \4gtgtt\d£di 

We estimate I 3 in exactly the same way as I2, which gives 

\I' 3 \ < A 3 \\gt\\ L 4gtt\\ 3 L2 
where A' 3 = 64M/3. We estimate I 3 in a similar way to I\ as 



< 12 sup 



which by use of (2.25) and the Cauchy-Schwartz inequality gives 



m<A!i\\ 9i \\ L A\9tt\\h 



where A' 3 ' = 64/3. 



Combining these estimates, we get (2.20) with 

A = 48+ M 

3 



(2.27) 



(2.28) 
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where M is the maximal-function constant in (2.27). □ 



Using ( |2.20[ ) in ( |2.18| ), we find that 

^ll.9«IU 2 < \ (2A \\9(,\\LA\9ii\\h- 
Since ||^||| 2 < ||g|| L 2 |lff«IU 2 and llflIU 2 = IMU* is conserved, we get 

^ I II s 1 2 ^11 |il/2|i ,,5/2 



(2.29) 



provided that (2.3) holds. It follows from (2.29) and GronwalPs inequality that if 



N < £o j where eo is sufficiently small, then ||<?f£||z,2 remains finite and ( |2.3[ ) holds in 
some time-interval < t < fc/e 2 , where the constants eo,k > may be chosen to 
depend only on ||ito||#- The same estimates hold backward in time, so this completes 
the proof of Theorem |1.1| 



By solving the differential inequality ( 2.29 1 subject to the constraint ( 2.3 1, we can 



obtain explicit expressions for eo and k. Let 

E = ||5olli2 4 ||5o«lli2 4 , 
which is comparable to ||uo||h 2 from (2.6). Then we find that Theorem 1 1 . 1 1 holds with 

k = 



eo = 



1 1 

2V2N E ' 



2 1 

3AE% 



where N is the constant in i2.2l and A is the constant in (2.28) 



3. Normal form transformation. In this section, we relate the near-identity 
transformation of the independent variables used above to a more standard normal 
form transformation of the dependent variables, of the form introduced by Shatah [7] 

v = u + B(u, u) 

where B is a bilinear form. 

We consider the normal form transformation given in [1] : 



le\d x \(h 2 ) 



H[u]. 



(3.1) 



Here, d x denotes the derivative with respect to x and \d x \ — Hd x . Differentiating 
(3.1) with respect to t, using ( |1.1| to eliminate u t , and simplifying the result, we 
find that this transformation removes the nonresonant term of the order e from the 
equation and gives 



v t + ^e 2 \d x \ [h\d x \(u 2 )] =H[«] 



(3.2) 



The bilinear form B in (3.1 1 is not bounded on H , but one can show that the normal 



form transformation (3.1) is invertible on a bounded set in H 2 when e is sufficiently 



small. We were not able, however, to obtain ff 2 -estimates for v from (3.2), because 



(3.2) contains second-order derivatives, rather than first-order derivatives as in (1.1) 



and there is a loss of derivatives in estimating the _ff s -norm of v. In fact, for every 
power of eu that one gains through a normal form transformation of the dependent 
variable, one introduces an additional derivative. 
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The appearance of additional derivatives is a consequence of using a zeroth-order 
linear term H[u] to remove a first-order quadratic term euu x . By contrast, higher- 
order linear terms lead to normal form transformations that are easier to analyze. For 
example, consider the KdV equation 

u t + euu x = u xxx . 

Then, assuming we can ignore difficulties associated with low wavenumbers {e.g. by 
considering spatially periodic solutions with zero mean), we find that the normal form 
transformation 



v = u 



leads to the equation 



G 



Vn 



In this case, the normal form transformation is bounded and it smooths the nonlinear 
term. 



To explain the connection between the normal form transformation (3.1 ) and the 



near-identity transformation (1.3), we reformulate (3.1| in a convenient way. Writing 



g = H[v] and taking the Hilbert transform of (3.1), we get the ODE 

g = h — ehh x . 



(3.3) 



We regard g(t,x) as a given function and use (3.3) to determine the corresponding 



function h. We may write (3.3) as 



h - 



hh x = 0, 



which agrees up to the order e with an evolution equation in e for h(t, x; e): 

h e -hh x = 0, h(t,x;0) = g(t,x). (3.4) 



By the method of characteristics, the solution of (3.4) is 
h[t,x;e) = g(t,£), 



which is the transformation (1.3). Since (1.3) agrees to the order e with a normal 



form transformation that removes the order e term from (|l.l|), this transformation 
must do so also, as we verified explicitly in Section [2j 



It is rather remarkable that the normal form transformation (3.1) can be implc 



merited by making a change of spatial coordinate in the equation for h, but we do not 
have a good explanation for why this should be possible. 
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